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Black hole accretion disks can form through the collapse of rotating massive stars. These disks 
produce large numbers of neutrinos and antineutrinos of electron flavor that can influence energetics 
and nucleosynthesis. Neutrinos are produced in sufficient numbers that, after they are emitted, they 
can undergo flavor transformation facilitated by the neutrino self interaction. We show that some 
of the neutrino flavor transformation phenomenology for accretion disks is similar to that of the 
supernova case, but also, we find the disk geometry lends itself to different transformation behaviors. 
These transformations strongly influence the nucleosynthetic outcome of disk winds. 
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I. INTRODUCTION 



Accretion disks are a compelling astrophysical setting, the properties of which are still being understood. Sophisti- 
cated numerical models have examined the origins of the disks and their structure. The disks may arise out of mergers 
between neutron stars , between a neutron star and a black hole 3] , or as a result of some stellar collapses [4H6| . 
Disks of the latter type emit primarily electron and anti-electron type neutrinos and relatively little mu and tau type. 
It has been suggested that accretion disks could play an important role in the production of gamma ray bursts M, |8| 
and they have also been studied as possible sites for nucleosynthesis in the absence of neutrino oscillations pjHl6j. 
Disks that originate from stellar collapse with sufficiently high accretion rates are understood to have a region of 
trapped neutrinos with mean energies of tens of MeV. The emission surface of the neutrinos generally exceeds that 
of the antineutrinos, but the antineutrinos have higher temperatures [Tol. [ItI [l8j. 

Beginning in the early 1990s, it was realized that coherent forward scattering of neutrinos could impact neutrino 
oscillation, resulting in coherence between neutrinos of different energies and parametric resonances [201 ] . Therefore, 
we expect that neutrino-neutrino interactions will play an important role in neutrino oscillations above disks. Neutrino 
oscillations involving neutrino-neutrino interactions have been extensively studied in general and in the contexts of 
the early universe and supernovae [2ll - [50l ] . When restricted to two flavors, the numerical results in regions of high 
neutrino density have been analyzed in terms of an analogy to precession and nutation of spins in a magnetic field [5ll - 
This analo gy i s called the Neutrino Flavor Isospin (NFIS) picture. Another approach was to use the techniques 
of BCS theory J47(. The early calculations assumed that interacting neutrinos shared a common history, in what 
is known as the single angle approximation [54j . Calculations that compute the histories of neutrinos across many 
trajectories resulted in effects not seen in the single angle approximation, including the decoherence of different energy 
modes (5f| [56[ . Numerical calculations [57l - [60l ] showed that oscillations split the neutrino spectra among flavors at 
discrete energies. These so-called spectral splits depend on the adiabaticity of the oscillations [13, IH|. They can be 
understood in terms of the NFIS picture as a magnetic resonance phenomenon (62|. 

Neutrino- neutrino interactions above disks were studied by Dasgupta et al. [63[. They generalized synchronized 
behavior among neutrinos of different energies to a disk geometry. For their disk geometry, Dasgupta et al. studied 
two identical, circular disks, without black holes at the center, emitting both neutrinos and antineutrinos, with 
number densities fixed to a constant ratio, so that neutrinos dominated the self-interaction term. Coherent forward 
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scattering of neutrinos on electrons (matter) was included by assuming a small mixing angle. With this model, a single 
"nutation" region, similar to that found in supernovae, was found when the neutrino densities became sufficiently low. 

In this paper, we create disk models that reflect the qualitative understanding we have of black hole accretion disks 
and include an explicit matter interaction term. We choose our neutrino and antineutrino disks to have different 
temperatures and different radii, so that the ratio between the neutrino flux densities will not be fixed, but will vary 
with position. 

Typically, the antineutrino disk will be hotter than the neutrino disk, but smaller. Therefore, there are two basic 
relationships the densities may have above the disks. Neutrinos may always dominate or antineutrinos dominate near 
the disk while neutrinos dominate further away. Our disks have "holes" in the center, from which no neutrinos or 
antineutrinos are emitted. This region corresponds to the space close to the black hole within its last stable orbit. To 
these disks, we add a nontrivial electron density. 

We also consider the role that neutrino flavor transformation plays in nucleosynthesis. The impact of matter- 
enhanced neutrino flavor transformation core collapse supernovae nucleosynthesis has been considered in e.g. [63467| 
and the impact of self interaction on supernova nucleosynthesis has been considered in [50l l68j . However, the impact 
of flavor transformation in disk nucleosynthesis has not been previously evaluated. 

In section [XT] we outline our two general disk models and in section IIIII we outline our method. In section IIV1 
we perform a single angle numerical computation of three flavor neutrino oscillations in these disks. In section |V1 
we analyze the results, and find a new type of transition region in accretion disks. In section IVII we discuss the 
ramifications for nucleosynthesis and in section fVIII we conclude. 

II. DISK MODEL 
A. General Features of Accretion Disks 

Our goal is to use a model that reproduces some general features of disks. The material near the plane of the disk 
is hot and protons and neutrons exist as free nucleons participating in the reactions of equation [T] as well as other 
scattering interactions. 

n + v P + p 

4-- -+. (1) 

For some disks with low accretion rates, below ~ O.f Mq/s, the reactions primarily proceed from right to left since 
electron neutrinos and electron antineutrinos are not trapped [69j |. 

As the accretion rate increases above ~ 0.1 Mq/s, the rates for neutrino interactions increase as well. Where 
the disk is hottest, near the plane of the disk, the neutrinos may be trapped [I?} • Further away, they can stream 
freely. For disks with high accretion rates, ~ 1 Mq/s, both electron neutrinos and antineutrinos are trapped [Tol.[T8|. 
although the emission surface of the electron neutrinos is more extended along the plane of the disk. In addition, 
electron antineutrinos begin to free stream closer to the plane of the disk than electron neutrinos. Therefore, electron 
antineutrinos are hotter than the electron neutrinos when they are both free streaming, but the electron neutrinos 
originate from a larger emission surface. 

Other flavors of neutrinos, including and v T are trapped in disks with high accretion rates, eg ~ 10 Mq/s 
[H, [lO] • They can produce a range of nucleosynthesis [7l[ and are interesting because they may be sites for r-process 
nucleosynthesis. We will not address them here, because they do not typically originate from stellar collapse. 

Given these considerations, we focus on two disks that describe two qualitatively different pictures. The flux above 
a disk can be antineutrino dominated close to the emission surface (since the antineutrinos are hotter) but neutrino 
dominated far from the disk since the neutrino emission surface is larger, and the disk is radiating net lepton number. 
This is the situation that one would expect from the type of steady state disk models described in Surman and 
McLaughlin [To| and Chen and Beloborodov in [ijj]. The other situation that can occur is that the flux above the 
disk is always neutrino dominated. 

B. Our Models 



We use two models that reproduce the features discussed in the section Hi Al We take disks to be thin, geometrically 
flat with a hole in the center with a constant temperature throughout. These temperatures determine the neutrino 
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Model 


T v 


Tp 


Rv 


Rv 


Ro 


A 


3.2 MeV 


3.4 MeV 


1.5 x 10 7 cm 


10 7 cm 


3.2 x 10 6 cm 


B 


3.2 MeV 


4.1 MeV 


1.5 x 10 7 cm 


10 7 cm 


3.2 x 10 6 cm 



TABLE I: Parameters corresponding to models in this paper, including T„, the temperature of the neutrino emitting surface; 
Tv, the temperature of the antineutrino emitting surface; R v , the outer radius of the neutrino emitting disk; Rv, the outer 
radius of the antineutrino emitting disk; and Rq, the inner radius of both the neutrino and antineutrino emitting disks. 

fluxes entirely since we characterize the spectrum with Fermi-Dirac distributions. The disks emit only electron 
neutrinos and antineutrinos and we assume that the neutrinos and antineutrinos emitted from the disks follow straight 
line paths, not relativistic paths. 

We chose two relationships between the neutrino and antineutrino disks. These disks are summarized in table HI 
First, we look at a disk where the neutrinos are slightly cooler than antineutrinos but have a larger disk. We call 
this model A, since it is the most qualitatively similar to the proto-neutron star: the space above this disks will 
always be dominated by neutrinos. We also look at a disk which is motivated by the disks discussed in [T3, EH where 
antineutrinos are significantly hotter than the neutrinos, but again the neutrinos have a larger disk. We call this 
model B, since the space above these disks will start out dominated by antineutrinos near the disk surface and at 
some further position will be dominated by neutrinos. 

We take the electron densities to be consistent with a wind type of outflow as discussed e.g. in [7lj |. and use the 
same electron density in both models. We give the explicit form of the parametrization in section IVT1 

We examine both disks in the normal and inverted hierarchies. The vacuum neutrino mixing parameters were taken 
to be the mass squared differences, to 2 — m\ — Sm^i = 7.59 x 10~ 5 eV 2 , |m 2 — m\\ = \Sm\ 2 \ = 2.43 x 10~ 3 eV 2 ; and 
the mixing angles, 6*13 = 9°, 6*12 = 34.4° and #32 = 45°. These choices are consistent with the current PDG values 

III. NEUTRINOS 

Neutrino oscillations are governed by an equation with the ambient neutrino potential, the vacuum contribution 
and a contribution from interactions with electrons and positrons, 

ij t S(n,q,x,t) = (H v (q) + ff e (x,t) + H w (x,t)) S(fi,g,x,t) (2) 
where t is the time elapsed, x is the neutrino position, and 

H VI/ (-x.,t) — V2Gf / dq (1 — cos6 M >) (p v (Cl, q,*.,t)dn v (Cl, g,x,i) — pa(Cl, g,x,i)dn s (f2, g,x,i)) (3) 
Jo Jn 

The integration variable, q is the energy of an ambient neutrino that shapes the potential. The infinitesimals, 
dn v (£l,q,~x.,t) and dna (O, q, x, t) are the densities of neutrinos and antineutrinos and <d qq ' is the angle between a 
neutrino with momentum q' and an ambient neutrino with momentum q. The information about density evolution is 
contained in S(Cl, q, x, t), the scattering matrix that evolves the neutrino density matrix, p v (Q, q, x, t): p v (£l, q, x, t) — 
S(Cl,q ) x,t)p I/ (p,,q ) x,0)S'(Q ) q,x. ) t). The antineutrino density matrix, Pu(yi, q, x, t) is similarly evolved by another 
scattering matrix, S(£l,q,x,t). We will take both the neutrino and antineutrino density matrices to start with all of 
the neutrinos and antineutrinos in the electron flavor. The matrix, Hy is the vacuum Hamiltonian. 

The matter Hamiltonian is proportional to T4(x, t) = \/2Gp{N e - (x, t) — N e + (x, t)), the contribution from neutrinos 
interacting with electrons and positrons, where N e ^ (x, t) is the electron(positron) density which is then scaled by 
V^Gf, where Gf is Fermi's constant. 

The difference between the accretion disk case and the protoneutron star is the geometry of the source. To handle 
the symmetries of the disk, we define a coordinate system based on the location of a neutrino, as shown in Figs. 
Q] and O We consider a neutrino that is a distance x from the center of the disk and a distance z above. From 
these coordinates, we define 9, the angle between the z axis and the neutrino momentum. We also define <j>, the 
angle between the x axis and the projection of the neutrino momentum. This information completely characterizes 
straight-line trajectories. 
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FIG. 1: A side view of the disk. The solid (red) star indicates the origin of the neutrino on a fiat disk with a straight-line 
trajectory. The empty star indicates the current position of the neutrino above the disk. 
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FIG. 2: A bird's eye view of the disk. The z axis is coming out of the page. The solid (red) star indicates the origin of the 
neutrino on the disk. 



In the coordinate system of Figs. [T]and[21 the angle between the ambient neutrino with coordinates <fi and 9, and 
the other neutrino with coordinates <p' and 9' can be written 

cos Q qq ' — cos 9' cos 9 + sin 9' sin 9 cos(0' — <j>) , (4) 

for all pairs of neutrino trajectories. The electron neutrino and antineutrino densities are G?n„(f2, q, x, t) and 
dn.£i(fi, q, x, t) respectively. 



dn„(p)(Q.,q,x,t) = — , (5) 

where is the flux of neutrinos(antineutrinos). 

In this work, we use the single angle approximation. This assumption factorizes the angular and energy integrals 
in equation [2] and is equivalent to assuming all neutrinos along every trajectory behave identically. Because we apply 
a single angle approximation and because we take the disks to have constant temperature, we can factor the integral 
of equation [5] into more manageable parts, 



dq (1 - cos9 gg -) (p u (Q, q,x.,t)dn v (Q, q,x,t) - py(fi,g,x,£)tZn*(f2,g,x,£)) 



y/2G F 

=-j- / ®v(q)dq / (l-cose M /)p„(f2,Q',x,t)#d(cos0) (6) 
47r Jo Jn v 

— -— / <5>v{q)dq ( (1 - cos Q qq >) p p (il, q, x, t)d(f>d(cos 9). 
47r Jo Jn B 
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where the fi„(^) reminds us to take the integral over the appropriate range of angles for the neutrino (antineutrino) 
disk. This term can be further simplified and the 9 integrals can be performed analytically. 

— G z f°° ( r R " f R " \ 

ff w (x, t) / *„(E)p u (E, x, t) / C(r, x, t)rdr - ® 9 (E)pl{E, x, t) / C(r, x, i)rdr dE, (7) 

V27T Jo y JRo JRo J 

where 

7r((l + m) (z cos 9' + x sin 0' cos 6') - Ax sin 0' cos 0' (a; 2 + z 2 )) 2E ( ^) 
C(r, x, t) = -2 ^ — ^ ^ U- , 8 

with i?o is the last stable orbit, R v (ja is the radius of the (anti)neutrino disk, E is the (anti)neutrino energy I — 
(x ~ r) 2 + z 2 , m = (x + r) 2 + z 2 and E(M) is the complete elliptic integral of the second kind. Therefore, we perform 
both the energy and r integrations numerically. 



IV. CALCULATION 



Using the disk models described in section |TT] and the technique described in section ITTll we calculate neutrino flavor 
transformation in both hierarchies. We follow the path of neutrinos from a position of x — z — 3 x 10 6 cm along 
the x = z trajectory. We plot the results in Figs. [3j [4j [5] and [6l The upper plots show the neutrino (red, solid) 
and antineutrino (red, dot-double-dashed) weighted survival probability as a function of progress along the trajectory. 
The weighted survival probability is taken to be the ratio between the oscillated and unoscillated neutrino capture 
rates, i.e. the rates of Eq. [TJ and is independent of energy. The lower plots show the interaction strengths for various 
contributions to the Hamiltonian in units of ergs. The red line is Tr{H ulJ ). If the term is dominated by neutrinos 
Tr(H uv ) > 0, the line is solid and if dominated by antineutrinos Tr(H vv ) < 0, it is dot-double-dashed. The cyan 
dotted line shows V e (t) while the blue dashed and green dot-dashed show the vacuum strengths, 

Smf 3 Sm 2 12 

2{E) y 2{E) 1 ' 

for the mean energy, (E) . 

First we use model A where the neutrino disk is larger than the antineutrino disk and slightly cooler, so that the 
neutrino flux is always larger than the antineutrino flux. We show the results for each hierarchy in Figs. |3]and|3J 
In the inverted hierarchy, dramatic transitions are first seen when electron density drops enough to come close to the 
Sm 2 3 vacuum scale just before 10 8 cm; in the normal hierarchy, the probability remains almost flat until transitions 
occur at ~ 2 x 10 8 cm. 

Our second model, model B, entails a neutrino disk that is large, but cooler than the antineutrino disk. The 
antineutrino disk contributes more to the flux close to its surface but further away, the neutrino disk takes over. 
The contribution to the evolution from neutrino-neutrino interaction changes sign somewhere above the disk. The 
results are shown in Figs. [5] and H3 The results are similar for both hierarchies. The weighted electron neutrino 
survival probability begins as unity. As the neutrino-neutrino interaction strength changes sign, the survival proba- 
bility plummets. The transition begins when the neutrino-neutrino interaction strength is approximately the same 
magnitude as the neutrino-matter interaction strength (up to the vacuum contribution) but opposite in sign. The 
survival probability rises for the antineutrinos again as neutrino-neutrino interaction strength and the neutrino-matter 
interaction strength are again the same. 



V. ANALYSIS 



We see three kinds of regions among our models, and associate them with different circumstances. 

(I) Standard MSW transitions that occur when the neutrino-matter interaction strength cancels the vacuum 
strength. 

(II) Transitions that occur when the neutrino-neutrino interaction strength is about the same scale as the vacuum 
strength. These are associated with nutation oscillations in the NFIS picture. 
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FIG. 3: Model A; inverted hierarchy. Weighted electron neutrino survival probability and strengths of the neutrino- 
neutrino interaction term for the inverted hierarchy. Upper plot: Weighted survival probability of electron neutrinos (red, 
solid) and electron antineutrinos (red, dot-double-dashed) as a function of the progress along the trajectory. The weighted 
survival probabilities are ratios of the energy integrated neutrino and antineutrino capture rates. All neutrinos began in the 
electron flavor at x — z = 3 x 10 6 cm. Lower plot: Comparison between the strength of the neutrino-neutrino interaction (red 
line), the vacuum strengths (blue and green dashed lines), and the neutrino-matter interaction strength (cyan dotted line) The 
neutrino and antineutrino terms never exchange prominence and the interaction strength is always positive. Regions of types 
(UJ) and (UU) are labeled in the upper panel. 



(Ill) Transitions that occur when the neutrino-neutrino interaction strength is the same size as the neutrino-matter 
interaction strength, and both are much larger than the vacuum strength. These are associated with the 
canceling of those terms. 

Standard MSW transitions, type (jj) are well understood. They become important to both our models beginning 
at several times 10 s cm in Figs. EHHEJ and El when low energy contributions from the neutrino flux begin to hit the 
MSW resonances. Neutrinos will undergo resonant transitions when the neutrino-matter interaction strength reaches 
the vacuum strength, 8ra\ 2 /2E; The hierarchy determines whether neutrinos or antineutrinos undergo oscillation 
enhancement at 8m \J2E. The MSW transitions are energy dependent, but since we have presented survival prob- 
abilities with the energy dependence integrated out, there is no single position where the effect takes place. In Fig. 
[3] for model A, a type (p} region associated with 8m\ z /2E occurs between about 2 x 10 8 cm and 10 9 cm, where the 
survival probability for electron antineutrinos rises. In the same figure, the region associated with 5m\ 2 /2E appears 
between about 10 9 cm and 10 10 cm where the electron neutrino survival probability drops slowly. Similarly, the type 
(HI region associated with Smf 3 /2E for model B in the inverted hierarchy is where the survival probability for electron 
antineutrinos drops between about 2 x 10 s cm and 10 9 cm in Fig. [5] In the normal hierarchy, the MSW effects are 
most important for 5m\ 2 /2E, again, appearing between about 10 9 cm and 10 10 cm. Figure U] shows a small rise in the 
electron neutrino survival probability in this region for model A, and Fig. [6] shows a significant rise in the electron 
neutrino survival probability for model B. 

A region of type (jlTj) appears in our model A in the inverted hierarchy just before 10 8 and in the normal hierarchy 
somewhat after 10 8 cm as seen in Figs. [3] and [4] The transformations occur as the neutrino-neutrino interaction 
strength approaches the vacuum scales, 8m\ 3 /(2(E)) and 8m\ 2 /(2 (El)). This causes a drop in probability for both 
neutrinos and antineutrinos in both hierarchies. Nutation oscillations in an anisotropic environment have been studied 
previously for the protoneutron star (UH^]. In the NFIS picture, as articulated in [53| and [3l|, strong oscillations 
are expected in such a region and are associated with nutation in the inverted hierarchy, where as smaller oscillations 
are expected in the normal hierarchy. 



7 




lxlO 5 lxlO 6 lxlO 7 lxlO 8 lxlO 9 



Progress (cm) 



FIG. 4: Model A; normal hierarchy. Weighted electron neutrino survival probability and strengths of the neutrino-neutrino 
interaction term for the normal hierarchy. Upper plot: Weighted survival probability of electron neutrinos (red, solid) and 
electron antineutrinos (red, dot-double-dashed) as a function of the progress along the trajectory. The weighted survival 
probabilities are ratios of the energy integrated neutrino and antineutrino capture rates. All neutrinos began in the electron 
flavor at x — z = 3 x 10 6 cm. Lower plot: Comparison between the strength of the neutrino-neutrino interaction (red line), 
the vacuum strengths (blue and green dashed lines), and the neutrino-matter interaction strength (cyan dotted line). Regions 
of types ([I]) and ([II]) are labeled in the upper panel. 



Previous calculations for neutrino flavor transformation in disks in the inverted hierarchy [63| followed along the same 
lines: both neutrinos and antineutrinos experienced nearly the same drop in probability in the inverted hierarchy in the 
type ((TTJ) region. We verify such a drop, at around 10 s cm in Fig. [3]but we find that neutrinos and antineutrinos behave 
differently. This is expected since we consider disks with neutrino and antineutrino trapped regions of different sizes 
(also with black holes at the center). Thus in our calculations, the effective angles for the neutrinos and antineutrinos 
are different. 

The normal hierarchy has not previously been considered with a disk geometry. However, it has been studied in a 
spherical geometry i.e. the protoneutron star. Bipolar oscillations were studied for the normal hierarchy by Duan, 
Fuller and Qian [55[ , who found small amplitude oscillations in regions of type (pi]) . The amplitude of these oscillations 
was found to depend on the effective mixing angle. The effective mixing angle in our type (|LT|) region of Fig. U is 
between 0.6 and 1 radians. This angle is much larger than those considered in [55j | and a larger angle is consistent 
with larger amplitude oscillations. 

In model B, we find little evidence for transitions that correspond with the type (JTTJ) region. 

Oscillations arise in a region of type as well. These occur in the normal and inverted hierarchies in model B, 
but not in model A. In model B, the region of type occurs at about 4 x 10 6 cm and then at 5 x 10 7 cm. To 
understand why this type of region can produce oscillations, we look at two flavors, simplifying the neutrino-neutrino 
interaction term to show the structure, 

ffw(t)= /X(t)-M*) h (t) y (10) 

The contributions from electron neutrinos, muon neutrinos, electron antineutrinos and muon antineutrinos are given 
by h„ e , h v , h i , c , hp respectively, so that eg, 



h Ve (t) = V2G F dq (1 - cos Q qq >) p u<ee (ft, q, x, 
Jo Jn 



t)<Mn,9,x,t), (ii) 
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FIG. 5: Model B; inverted hierarchy. Weighted electron neutrino survival probability and strengths of the neutrino- 
neutrino interaction term the inverted hierarchy. Upper plot: Weighted survival probability of electron neutrinos (red, solid) 
and electron antineutrinos (red, dot-double-dashed) as a function of the progress along the trajectory. The weighted survival 
probabilities are ratios of the energy integrated neutrino and antineutrino capture rates. All neutrinos (antineutrinos) began in 
the electron flavor at x = z — 3 x 10 6 cm. Lower plot: Comparison between the strength of the neutrino-neutrino interaction 
(red line), the vacuum strengths (blue and green dashed lines), and the neutrino- matter interaction strength (cyan dotted line). 
Regions of types ((I]) and are labeled in the upper panel. 



where we've suppressed the dependence of h(i) on x. They correspond to the probability of finding a neutrino or 
antineutrino in a particular flavor state up to the rest of the interaction scale. For example, at the start of our 
trajectory, before any oscillation, there are only electron flavor neutrinos and antineutrinos, so that h Vfl = h Pfi = 
and 



h v ,(t) - h B ,(t) = V2Gf ! dq f (l-cosG, 
Jo Jo. 



>) (dn„(fi, q, x, t) - dn p (Q, q, x, t)) . (12) 



At an arbitrary time, the total Hamiltonian will be 
H v + H e (t) + H vu (t) = 

1 / Ve (t) h »-(t)-^(t) + ^.W-^W _ cos2Qi2 sg. sin26i2 + 2hefi(t) \ 

2 ^ sin2012 + 2 Ke {t) eg cos 26 12 - V e (t) + "- ( "-^ ( " - feBe(t) ; feg " (t) J ' 

(13) 

where we've subtracted of the trace of the interaction terms which supplies an overall phase. Like with the MSW 
resonance, we expect to see large flavor transformation when the diagonal terms vanish, leaving only the flavor mixing 
off-diagonal terms remaining. When there are only electron flavor neutrinos and antineutrinos, the interaction term 
is 

H v + H e (t) + H uu (t) = 

i(v e (t) + ^--^-eg-cos29 12 eg sin2 e 12 + 2h eil (t) \ (u) 

2 ^ s™. sin 2t?i2 + 2h tie {t) cos2c? 12 - V e (t) - ^ + ^ J ' 

Neutrinos and antineutrinos will experience the cancellation slightly differently because of the relative sign of the 
vacuum portion to the interaction portion. Nevertheless, since the neutrino-matter interaction strength is very large 
compared to the vacuum, there will be transitions for both when the system is antineutrino dominated, h„ e (t) > h Vt , (t). 
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FIG. 6: Model B; normal hierarchy. Weighted electron neutrino survival probability and strengths of the neutrino-neutrino 
interaction term for the normal hierarchy. Upper plot: Weighted survival probability of electron neutrinos (red, solid) and 
electron antineutrinos (red, dot-double-dashed) as a function of the progress along the trajectory. The weighted survival 
probabilities are ratios of the energy integrated neutrino and antineutrino capture rates. All neutrinos (antineutrinos) began in 
the electron flavor at x — z = 3 x 10 cm. Lower plot: Comparison between the strength of the neutrino-neutrino interaction 
(red line), the vacuum strengths (blue and green dashed lines), and the neutrino-matter interaction strength (cyan dotted line). 
Regions of types ^ and (|III|) are labeled in the upper panel. 

These criteria are met in model B at about 4 x 10 6 cm. We see large flavor transitions for both neutrinos and 
antineutrinos. When there are no electron flavor neutrinos, the interaction term becomes 

H v + H e {t)+H vv {t) = 

1 / Fe ( t )_^W _4^ COS 20 12 ^Bin2fl ia + 2M*) ^ (15) 

2^ 4^sin2# 12 + 2VW ^co S 29 12 -V e (t) + ^-^ J ' 

In this case, we expect that there will be big transitions in both hierarchies when the system is neutrino dominated, 
h v (t) > hp (t). These criteria are met in model B (Figs. [5]and[H]) at about 4 x 10 6 and 5 x 10 7 cm and we do indeed 
see flavor transitions, this time, larger for antineutrinos and smaller for neutrinos. Note that we do not see large 
flavor transitions in model A at about 10 6 cm because the neutrinos at that point are all electron flavor and neutrinos 
outnumber antineutrinos. There is no way for the diagonal terms to cancel in model A, until the MSW region at a 
few times 10 8 cm. 

In model B, a region of type (JTTJ) occurs during the first type (|III[) region. However, we associated this transition 
primarily with a type (|III[) region since the starting point corresponds to the approximate cancellation of the diagonal 
term in the Hamiltonian. 



VI. RAMIFICATIONS FOR NUCLEOSYNTHESIS 

The neutrinos emerging from the disk will interact with material outflowing from the inner disk, and thus influence 
any resulting nucleosynthesis. In this work, we present an example case to illustrate the potential impact. We begin 
with a one-dimensional disk model with rh = 3M Q /s, a — from [18j : this disk has trapped regions for the electron 
neutrinos and antineutrinos, but not the mu and tau flavors. We calculate the neutrino and antineutrino decoupling 
surfaces as in • We then find the single-temperature, flat disk approximations to these surfaces that produces the 
best match to the full disk neutrino emission. The resulting parameters — T„ = 3.2 MeV, R v = 150 km, T p = 4.1 
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MeV, Rp = 100 km — are identical to those of model B, where antineutrinos dominate near the disk and neutrinos 
farther out. 

The element synthesis is calculated as in [7l|. The outflow is taken to be adiabatic and radial, with velocity 
v = i>oo(l —To/t)^, where Voo is the final coasting velocity of 0.1c, and controls how rapidly the material accelerates. 
Close to the disk the the material consists of primarily neutrons and protons, with a ratio set by the forward and 
reverse weak reactions, Eq. [1] The subsequent nucleosynthesis is calculated with the sequence of network codes 
described in [7l[; neutrino interactions on nucleons are included throughout. 

If neutrino oscillations are ignored, the outflowing neutron-rich material can synthesize nuclei characteristic of the 
second abundance peak of the r-process of nucleosynthesis, as shown in Fig. [7] for outflow parameters s/ks = 50 
and /3 — 1.4. The initial neutron-to-proton ratio in the material is greater than one, since the antineutrino capture 
on protons is initially favored over neutrino capture on neutrons. All of the protons and most of the neutrons 
assemble into alphas and heavier 'seed' nuclei; the r-process proceeds when the remaining neutrons are captured on 
the seeds. However, once the alphas form, the remaining neutrons are depleted by neutrino interactions on neutrons, 
and additional alphas form immediately from the protons that result. Thus, there are too many seed nuclei and too 
few remaining neutrons for a full r-process. The situation changes when neutrino oscillations are taken into account. 
The disappearance of neutrinos after they have progressed 4 x 10 6 cm causes the neutrino interaction rate to plummet, 
as shown in the top panel of Fig. [5] Thus, the conversion of neutrons to protons during alpha particle formation is 
cut off. Fewer alpha particles and seed nuclei form, and more neutrons remain to capture on these seeds, as illustrated 
in the bottom panel of Fig. [5J The result, shown in Fig. [7J is a more robust r-process that produces heavy nuclei out 
to the A ~ 195 region. 

As the top panel of Fig. [8] shows, the neutrinos reemerge at around 1000 km, and thus could potentially slow the 
r-process by converting neutrons to protons at later times in the nucleosynthesis. To see the extent to which this 
effect operates, we additionally run the simulation with the neutrinos turned off after the neutrinos have proressed 
> 2 x 10 7 cm. The results of this simulation are included in Figs. [7Jand[5J which show that by this time the capture 
rates are small and so the impact on the r-process is minimal. 

VII. CONCLUSIONS AND REMARKS 

We have studied models of neutrino and antineutrino emission for accretion disks that encapsulate the qualitative 
behavior of the neutrino fluxes leaving the disk. Two models were examined in detail: those dominated by neutrinos 
(model A) and those that begin dominated by antineutrinos and end up dominated by neutrinos (model B). The neu- 
trino dominated disks in the normal hierarchy result in little flavor transition until the neutrino interaction strengths 
become close to the vacuum strengths. In both hierarchies they exhibit oscillations in type [TT] (nutation/bipolar) 
regions. On the other hand, disks that begin antineutrino dominated and end up neutrino dominated produce large 
flavor transition when the neutrinos flux is about the same as the antineutrino flux. These transitions are associated 
with the cancellation of the neutrino and antineutrino terms with the neutrino-matter interaction strength, i.e. a type 
IIHI (matter-neutrino enhanced) region. 

These transitions close to the disk are particularly important because they occur at a time when the neutrinos 
are influencing nucleosynthesis. We have shown the illustrative example of how neutrino oscillations in the normal 
hierarchy impact nucleosynthesis. Using a disk which approximates the type of disk found in a "collapsar" scenario 
Q, i.e. one that has trapped electron neutrinos and antineutrinos only, with different sized trapping regions, we find 
that the addition of neutrino oscillations enables the formation of r-process elements. This disk can produce r-process 
elements when neutrino oscillations are taken into account. 

Multiange effects have been shown to delay the onset of type II flavor transformation in core collapse supernovae 
(45j . and this delay can influence nucleosynthesis [SOj . Neutrino halo effects have been suggested as a mechanism to 
alter the simple picture of Type II transformations as well [74j]. Since some of the transitions discussed here occur 
close to the disk, multiangle oscillations should be considered in the future for their impact on type III regions as well. 
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FIG. 7: Final abundance patterns for simulations with disk outflow parameters s/ks = 50 and /3 = 1.4 in the absence of 
neutrino oscillations (red line), with oscillations assuming a normal hierarchy (blue line), and with oscillations until r > 2 x 10 7 
cm, when the neutrino interactions are turned off (green line). The scaled solar r-process pattern [73j is included for reference 
(black crosses). 
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